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Abstract 

We calculate the bispectrum of the gravitational field fluctuations generated 
during warm inflation, where dissipation of the vacuum potential during in- 
flation is the mechanism for structure formation. The bispectrum is non-zero 
because of the self-interaction of the scalar field. We compare the predictions 
with those of standard, or 'supercooled', inflationary models, and consider 
the detectability of these levels of non-Gaussianity in the bispectrum of the 
cosmic microwave background. We find that the levels of non-Gaussianity 
for warm and supercooled inflation are comparable, and over-ridden by the 
contribution to the bispectrum due to other physical effects. We also con- 
clude that the resulting bispectrum values will be undetectable in the cosmic 
microwave background for both the MAP and Planck Surveyor satellites. 

PACS number(s): 98.80.Cq, 98.70.Vc, 98.80.Es 

I. INTRODUCTION 

The most commonly adopted model of the early Universe pictures the large scale struc- 
ture present in the Universe today to be seeded by small scale fluctuations in the matter 
distribution of the primordial universe during an inflationary phase. In order to realize 
an inflationary regime, the generic dynamical model is based on a single scalar field, often 
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termed the inflaton. The homogeneous, zero-mode component of this field is pictured to roll 
on a ultrafiat potential, thereby sustaining a large potential energy and negligible kinetic 
energy, which are the necessary conditions for realizing inflation. These models assume there 
is no radiation production during the inflation period, with any radiation present prior to 
inflation rapidly diluting away. We refer to this picture as supercooled inflation, describing 
specifically the thermodynamic state of the universe during inflation. In this picture, the 
supercooled inflation phase is separated from the radiation-dominated regime by a brief 
reheating period, in which vast amounts of radiation are rapidly produced. 

In supercooled inflation the initial seeds of density perturbations result from quantum 
fluctuations of the inflaton fleld. To a good approximation, the fluctuations have a Gaussian 
distribution and produce a nearly scale-invariant spectrum. A perfect Gaussian distribution 
would imply that the density perturbations have no connected correlations higher than the 
2-point correlation function in real and Fourier space. However, the self-interaction of 
the inflaton fleld is known to produce non-zero, but extremely small, non-Gaussian effects, 
and there have been predictions calculated of these effects [|l| and their detectability in the 
distribution of the cosmic microwave background (CMB) fluctuations 0,0. Other structure 
formation scenarios, such as the class of defect models, and multiple-fleld inflation models, 
generally give larger deviations from Gaussianity. 

There are many ways of testing the Gaussian hypothesis, such as the genus and Euler- 
Poincare statistic 0-0], studies of tensor modes in the CMB excursion set properties 
PJTOd, peak statistics 



nHT^ and wavelet analyses (e.g. HTSHT^ 



Then there are the set 
of higher-order correlation functions, such as the three-point function (e.g 



UMI]), the 

bispectrum pT|-p^, and the trispectrum A signiflcantly non-Gaussian signal in the 



CMB sky, measured from the data of the Cosmic Background Explorer (COBE) Differential 
Microwave Radiometer (DMR) satellite instrument, launched in 1989, has been claimed P3 



but doubts have been cast on this as a signiflcant primordial signal, and more recent papers 
using this data [pi-|2H| and, on a smaller angular scale, the MAXIMA balloon experiment 
data of 1998 [0,^ flnd results largely consistent with Gaussianity. 

Gangui et al. ^ calculated predictions of the bispectrum for several variants of the 
supercooled inflationary scenario. In the limit of no instrument or sampling noise, the 
minimum variance on CMB data is the cosmic variance - arising from fact that our 
Universe could only be one of a Gaussian ensemble. Gangui et al. [|1| found the values for 
skewness for the bulk of their single fleld models to be considerably smaller in magnitude than 
that resulting from cosmic variance. The skewness is, however, only a single statistic which 
is part of a much wider class - the three-point correlation function. There is the prospect of 
doing much better by using the whole class or, equivalently, its harmonic counterpart, the 
bispectrum. 

COBE had an angular resolution of 7°. For COBE it has been possible to analyse all the 
modes of the bispectrum up to the resolution limit [^,^. The Planck Surveyor satellite is 
planned for launch in 2008. It will have an angular resolution of down to 5 arcminutes. The 
Microwave Anisotropy Probe (MAP) is operating now, with a resolution of 12.6'. Both of 
these experiments allow vast numbers of bispectrum modes to be analysed in principle, and 
it is interesting to see if either could distinguish warm inflation from supercooled inflation 
on the basis of the bispectrum. 

In this paper we predict the form of the fluctuations and quantify the non-Gaussianity, 
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using the bispectrum as a measure, for warm inflation dynamics. Warm inflation differs 
from the standard, supercooled picture of inflationary cosmology in that the process of ra- 
diation production becomes an important constituent of the theory. In particular, radiation 
production occurs concurrently with inflationary expansion and the presence of this radition 
influences the seeds of density perturbations. 

The warm inflation picture of inflationary dynamics is a comprehensive set of possible 
interactions between fields during inflation. In this picture no a priori assumptions are made 
about multi-field interactions, thus particle production, during the inflationary epoch. As 
such, the warm inflation picture makes explicit that the thermodynamic state of the universe 
during inflation is a dynamical question. In particular, the supercooled inflation emerges as 
one limiting case in which the interactions are negligible. More commonly dissipation is pos- 



sible in warm inflation with a resulting density of radiation present during inflation pl|-|33 



Here both strong |l3ll^3| and weak [^,0,^ dissipative regimes have been examined which 
offer several variants to the basic picture. 

A variety of warm inflation models have been developed at a phenomenological level 
33| , |3^j36| , |37[] . From this it has generally been understood that warm inflation can solve 



the basic cosmological puzzles of horizon, flatness and density fluctuations. However, up to 
now no study has been made on the degree of non-Gaussianity typically emerging in warm 
inflation models and it is important to quantify such effects. 

In this paper, a general methodology is developed for computing non-Gaussian effects 
in warm inflation scalar field models. The statistic we use to quantify the resulting non- 
Gaussianity in the cosmic microwave background predicted for the case of warm inflation 
models is the bispectrum. This formalism is then applied to the A0^, and m^0^ models 
for the strong dissipative regime of warm inflation in the interest of discovering whether 
these effects are comparable, or indeed distinguishable, from the predictions of supercooled 
inflation. 



II. WARM INFLATION DYNAMICS 

The equation of motion for the zero mode of the scalar inflaton field, 0, in general has 
the form 

4> + 3H^ + T,p + V\^) =0, (1) 

where the overdots represent time derivatives, H = ^ is the Hubble parameter, and R{t) 
is the cosmic expansion factor. The necessary condition for inflation is domination of the 
inflaton potential energy over all other energy components in the universe. This is achieved 
by requiring the inflaton potential to have sizable magnitude and be very flat. The flat- 
ness of the potential allows slow roll motion of 0, so that inflaton kinetic energy becomes 
neglible with the effect that the term can be dropped from the equation of motion. In 
supercooled inflation dissipative effects, which in Eq.(ll|) are symbolically expressed through 
the term r0, are assumed negligible during the inflation period and only emerge during the 
subsequent reheating period. On the other hand, the basic observation of warm inflation is 
that inflationary conditions remain energetically possible even in the presence of a sizable 
radiation component. Thus in warm inflation dissipation remains active during inflation. 



3 



with the simplest representation of these effects being the form in Eq.(p. The resulting 
evolution equation for the inflaton in warm inflation therefore is 

d(j) 1 dV{(f)) 



dt 3H + T 



(2) 



where similar to supercooled inflation, the slow- roll condition is required, (3if + r)|(/)| ^ |(/)|. 

The presence of dissipation implies throughout the inflation period radiation is produced 
from conversion of vacuum energy. In the context of Friedmann cosmology, the stress energy 
conservation equation consisting of a radiation, pr, and a vacuum energy component, p„, is 

Prit) = -4:Prit)H ~ p,{t). (3) 

If there were there no dissipation, then p^ = and the radiation component would be rapidly 
red-shifted away as pr ~ e~^^*. However with dissipation, radiation is being produced 
continuously from conversion of scalar field vacuum energy. For dissipation of the form in 
Eq.®, -p, = r02. 

Note that in general the conversion of vacuum energy into radiation will result in some 
type of reaction back upon the scalar field, with the term being the simplest phenomeno- 
logical possibility. Attempts to obtain first principles scalar field evolution equations for 
warm inflation have obtained dissipative effects which in general are temporally nonlocal 
although limits also have been obtained in which the effects are of the local form in Eq.(|l]) 



31|,|38|J 

The production of radiation during inflation in general will influence the seeds of density 
fluctuations. In particular, if the temperature during inflation is bigger than the Hubble 
parameter, T > H, then structure formation can be significantly affected by the thermal 
component. In a preliminary work leading to the development of the warm inflation scenario 
||32|| , it was shown that a tiny dissipative component F> 10~^if already is adequate to realize 
T > H. Subsequently the phenomenology of both weak T < H [p2| , ^ and strong T > H 
3|,|38|-|0l dissipative regimes have been studied for warm inflation. More attention has 



been given to the strong dissipative regime. This primarily is because the main focus of 
first principles quantum field theory studies of warm inflation [|38| , ^P1 , ^^2| -|4^ have been 



in this regime, since it is the more difficult of the two, and once this regime is understood 
the weak dissipative regime easily would follow. In this paper, expressions will be obtained 
which apply to the strong dissipative regime, such that Eq.(^ reduces to ^ = —jr, and 
during inflation pr ~ 0, iprli ~ — p„. 

In order to treat the fluctuations of the inflaton field 50(x, t), it is assumed that they 
are small and the full inflaton field is expressed as 0(x, t) = (j)o{t) + 50(x, t) with 0o being 
the homogeneous 'background' field, (50(x, t) ^ (t>o(t)- The equation of motion for the 
fluctuations of the inflaton field can be obtained by imposing a near-thermal-equilibrium, 
Markovian approximation, which therefore implies the fluctuation-dissipation theorem is 
applicable. From this the equation of motion for the full inflaton field emerges as 



d(f){x,t) 1 



e 



-2m 



V^<j){x,t)-V'{(j){x,t))+7]{x,t) . (4) 



dt F 

Implementing the fluctuation-dissipation theorem immediately determines the properties of 
the noise. With respect to physical coordinates and in momentum space, these properties 
are 



4 



(^) = o 



(5) 



{rjik, t)7]{k', t')) = 2TT{27rf6^^\k - k')6{t - t') (6) 

III. THE STATISTICS OF WARM INFLATION PERTURBATIONS 
A. The Predictions and Properties of Gaussian Fields 

Single field inflation models broadly predict Gaussian primordial density fluctuations. 
Multiple field inflation models may lead to a non-Gaussian distribution (e.g. chi-squared) . 
When second-order effects are taken into account, however, there are corrections to these 
general predictions. There is a resulting non-Gaussian signal in the CMB, the magnitude 
of which varies depending upon the self-interaction of the inflaton field |ll9| , |46| ,p]| . 

The statistical properties of a Gaussian field with mean zero are fully contained in its 
power spectrum or its two-point correlation function in real space. Higher order correlations 
can be expanded in terms of these quantities, so correlation functions of even order of the 
distribution can be written as products of two-point functions, and correlations of odd-order 
can be written in terms of products of two-point functions and the expectation value of the 
field. As there are no connected correlations over the two-point, all odd order correlation 
functions of a multivariate Gaussian distribution with zero mean are equal to zero. For a 
non-Gaussian field, the higher-order connected correlation function can be non-zero. 

The three-point correlation function of the density perturbation distribution in Fourier 
space, otherwise known as the bispectrum pTl-p^,p]1 is the quantity we have chosen to evaluate 



as our measure of the non-Gaussianity generated by warm inflation. This quantity translates 
to the harmonic bispectrum of the CMB. 

B. The Warm Inflation Bispectrum 

We will begin by expanding Eq.(^) in order to obtain the evolution equations up to 
second order in the fluctuations (50(x, t) = (50i(x, t) + (502(x, t), where 6(f)i = 0{6(j)) and 
(502 = 0{6(j)'^). All calculations that follow will be in momentum space, and with respect to 
the physical versus comoving momenta, where recall 

Hereafter for physical momenta, the notation will have no subscripts fcphys = k, and mag- 
nitudes will be denoted without boldfacing as k = \k\. The choice of physical coordinates 
arises since we are interested in the evolution of the inflaton mode while they are sub-horizon 
scale. During this time, the dominant effect on the modes is from the high temperature heat 
bath, and these effects are more conveniently analyzed in physical coordinates. When evalu- 
ating the equation of motion for the fluctuations, the time dependence of the physical modes 
will be treated adiabatically with respect to the characteristic macroscopic time scale, the 
Hubble time ~ Thus the evolution equations to be written for the inflaton modes will 
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only be valid over a time interval ~ 1/H and a complete solution over longer time for a 
given mode can be obtained by piecewise construction. Thus the equations of motion of the 
first and second order fluctuations over a time period ~ are 



+V{k,t)] 



(8) 



dt 



r 



W"'{utmiik,tr] 



(9) 



where 6(f){k,t) = 6(j)i{k,t) + 6(j)2{k,t) is the infiaton mode with physical momentum k at 
cosmological - corresponding to the homogeneous background - time t. The properties 
of the noise are given in Eqs.(^|D. We also assume here that the evolution of scalar field 
fluctuations can be studied in a particular gauge where metric perturbations can be neglected 
compared with those of the scalar field itself. The latter assumption, originally made in Ref. 
|T9| , allows one to focus on the computation of those non-Gaussian features which are 
directly produced by non-linearities (i.e. self-interactions) of the scalar field itself, rather 
than by their backreaction on the underlying geometry. As we will see, this approach will 
give rise to a level of non-Gaussianity which is comparable to that obtained in Ref. [|l|, 
where backreaction effects were simply modelled through a local modification of the Hubble 
expansion rate during inflation. 

Dividing cosmic time into successive time intervals of order 1/H, tn — = 1/H, the 
solutions of Eqs. (H) and (|) for < t < tn are respectively 



5Mk,t) 



,v{k.t') 



A{k, t - tn^i) I dt''-^^^A{k, t' - t„_i)-i 

Jtn-l i 

+A{k, t - t„_i)50i(fce-^(*"-*"-), t„_i) 



(10) 



where 



Hk,t) = 
A{k,t-tn^i) f dt' 
dp^ 



B{t') 



6MP,t')SMk - P,t') 



(27r)3 

+A{k, t - t„_i)502(fce-^(*"-*"-\ t„_i) 



A{k,t' -tn-i)-^ 



A{k, t) = exp 



(11) 



(12) 



Bit) 



(13) 
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In both solutions Eqs. ([T0| ) and (pi]), the second term on the RHS are "memory" terms 
that reflect the state of the given mode at the beginning of the given time interval. The 
relevance of these memory terms leads to the important concept of freeze-out [^. By 
definition of freeze-out, for \k\ ^ kp the memory terms damp away within a Hubble time 
and for \k\ ~ kp they do not. To quantify this criterion, the freeze-out momentum kp is 
defined by the condition 

In general for warm inflation V"{(j)o) < TH, so the above condition can be simplified to 
kp = yVH. In supercooled inflation the freeze-out wavenumber would correspond to the 
Hubble scale, as the quantum fluctuation becomes classical on horizon exit. For warm 
inflation the fluctuations freeze in before horizon exit. 

Freeze-out implies noteworthy features about the solutions. When k > kp, since the 
memory terms are negligible, S(j)i, 502 primarily are determined by the state of the envi- 
ronment nearby in time. Then, when k < kp, 6(j)i, 5(j)2 are determined dominantly by their 
state at time of freeze-out. These two facts imply the time-slicing approach we used to 
solve for 501, 502 is well justified. 

From the solutions Eqs. (|T0|) and (0), we are interested in computing the three-point 
correlation function of the inflaton fluctuations at the largest observable scales, which in 
particular cross the horizon in the interval ~ 50 — 60 e-folds before the end of inflation. 
During this time interval, we will compute the defining parameters of the three-point cor- 
relation function, which as will be seen can be expressed in terms of a bispectrum that 
well approximates the Sachs-Wolfe regime, / ~ 50. Thus from the three-point function 
(50(fci, ^60)50(^2, ^60)50(^3, teo)), the amplitude and slope of the bispectrum are determined 
at time t ~ teo, 60 e-folds before the end of inflation and for ki, k2, k^ all all within a few 
e-folds of exiting the horizon. 

The leading order contribution to this three-point correlation function comes from two 
first-order and one second-order fluctuation as 



(50(fcl,t)50(fc2,t)50(fc3,t)) = 

A{k^, t - teo - l/H) / A-\k,, t' - t,o - l/H)B{t') 
n{ki, ti)50i(p, t')){6Mk2, ^2)501(^3 - P, t')) 



(2^3 
+A{k3,t-teo 

502(fc3e~\t6O ■ 

^ ^3)) ^ 



-l/i/)(50l(fcl,t6o)50l(fe2,t6o) 

- l/H)) 
{k2 ^ ks)) 



(15) 



In this expression, since B{t') is slowly varying, it can be approximated as a constant. 
Similarly S(j){k, t) can be fixed at its freeze-out value. The three-point function on the RHS 
arises from the memory term of 502 in Eq. (|TID, since k < kp. In evaluating this quantity, 
first note the coefficient in front can be approximated as unity, A{k,tn-i — tn-2) ~ 1- 
Furthermore, since all three momenta will overlap in the freeze-out region, the three-point 
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correlation function at t^Q — 1 / H is approximately the same as at t^o, and this property will 
repeat itself for the time interval 

AtF = tH-tp^j^ln{^), (16) 

where represents the time at Hubble crossing of the smallest of the three inflation pertur- 
bation modes, and tp represents the time when the last of the three wavevectors thermalizes. 
Thus Eq. ( |T5|) becomes 

{S(j){ki,t)d(j){k2,t)6<f){k3,t)) ^ B{teo)AtF 



+ (fcl ^ fcs)) + (fc2 ^ fcs)] . (17) 



C. Estimating the Magnitude of the Non— Gaussianity 

It is possible to write a general expression for the bispectrum for slow roll, single field, 
supercooled inflation models as well as for the set of warm inflation models. 

($(fci)$(fc2)$(fe3)) = 

Au,i{2Trf6^{k, + fcs + fcg) [P$(fci)P$(fc2) + perms] (18) 

The relation between the scalar field fluctuation and the gravitational field has the simple 
form [H7 



$(fc) = -^S-sm, (19) 

5 6 



thus y4inf for a strongly dissipative warm inflation regime is 

(20) 



.warm _ f <^ 

~ 3 [h 



H \H r 



Comparative estimates for the non-Gaussianity of the cosmic microwave background can 
then be calculated using A^^^ and the shape of the inflationary potential. To estimate the 
magnitude of AJ^f^™, consider the model 

v{ct>) = (21) 

in the region < < M, where, fixing the origin of time at t = 0, initially 0(0) = M. 
Here the self-coupling constant A is dimensionless. M sets a basic scale such as the Grand 
Unified scale ~ lO^^GeV, although the final answer for Af^^"^ is independent of this scale. 
The solution of the zero-mode evolution equation Eq. (H) for this case is 
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In [33 1 it is shown the number of e- folds of inflation for this quartic potential is 



1 / 487rr2^ 

The fluctuations in the scalar field caused by thermal interactions with the radiation field 



are 



40 



The temperature, T can be calculated using the relation = [g^n'^ /30)T^, where is the 
number of relativistic fields ~ 150, taking into account the relation between the radiation 
energy density and the scalar field potential presented in Eq. The CMB amplitude is 
given by 

Sh = IS-S<P, (25) 

o 



where from COBE data it is measured to be 6h = 1-94 x 10^^ p9|,|48|j49[ . 

Setting Ne = 60 and using Eqs. (p2|), (p3D and (^), the value of A is found to be 

A = 7.2 X 10^1^ (26) 



and this gives us all the quantities needed to evaluate Ajnf for warm inflation, Eq. (|20| ) 

AY^f^ = 7.44 X 10^1 (27) 

The equivalent quantity to Ai^f for supercooled inflation appears in Gangui et al. |^ as 
$3 = ^f'^P«^''™°i*^'i = 5_5g X 10~^ for a quartic potential. A similar quantity, f^L^ appears in 
30| and a related quantity in ||^, Ai^^. They are related as 

$3 = Ai„f = 2/jvL = A„fi(27r)6 . (28) 

For this potential we see that the non-Gaussianity in the curvature in warm inflation models 
is comparable to that in supercooled inflation. Furthermore due to the slow-roll behavior, 
for general models with local behavior of the form 0"? with 2 < g < 4, the value of A^^^'^ 
will be within the same order of magnitude as Eq. (pTj). Note, however, that this is only 
the contribution arising from non-Gaussianity in the inflaton field. In addition, there is a 
contribution at the last-scattering surface arising from second-order gravitational perturba- 
tion theory. Even at this relatively early time, the nonlinear gravitational equations induce 
a contribution Af^/^ "■'d^'- = C(l), [^§U[^J5ggg, which is therefore the dominant source 
of non-Gaussianity in both warm inflation and supercooled inflation models. The prospects 
of measuring this are poor; the bispectrum of the microwave background radiation is the 



most obvious possibility, and this has been investigated by |30]. Note that this assumes that 



the Hubble parameter is unchanged at horizon exit for all modes considered. The prospects 



are best for the Planck satellite [Q, as this measures modes up to £ ~ 2000, but even with 
the optimistic assumption that all foreground contaminants can be removed perfectly, the 
expected error on f^i is still 5. Since Silk damping causes the power to decline rapidly on 
scales smaller than the Planck beam, there is little prospect of a future, higher-resolution 
experiment being able to improve significantly on this. Similarly, polarisation measurements 
are unlikely to help, since most of the power is expected at £ ~ 100, for which there are 
relatively few bispectrum coefficients. 
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IV. CONCLUSIONS 



We have explored the evolutionary behaviour of gravitational field fiuctuations generated 
by warm infiation, up to second order in the scalar field perturbations, for the particular 
case of strong dissipation. We have obtained predictions for the non-zero bispectrum of the 
gravitational perturbations due to the self-interaction of the infiaton field, and the result- 
ing harmonic bispectrum for the Sachs- Wolfe region on the cosmic microwave background. 
Eq.(p!8D represents the gravitational bispectrum in a generic form, with the value of Ai^^ 
resulting from the theory of the fiuctuation generation mechanism. The CMB bispectrum 
can also be related to Ai^f, provided that the potential is sufficiently fiat that the Hubble 
parameter etc are approximately constant on horizon exit for all modes considered. We find 
that the inherently classical mechanism for the generation of fiuctuations in warm infiation, 
which is palpably different from the corresponding mechanism for supercooled inflation, pro- 
duce a level of non-Gaussianity of approximately the same magnitude, Ainf ~ 10~^. This 
arises from non-Gaussianity in the inflaton fleld itself. The dominant contribution to the 
curvature bispectrum, however, comes from general relativistic second-order perturbation 
theory, which contributes an Ai^f ~ 1, in both warm and supercooled inflation. Unfortu- 
nately this level of non-Gaussianity in the CMB bispectrum appears to be unobservable, 
even with the Planck satellite. Conversely, any measure of primordial non-Gaussianity sig- 
nificantly in excess of /at/, = 1 would rule out both warm infiation and supercooled infiation, 
although some alternatives such as the curvaton model [BB[ or multi-field infiation models 



e.g. ||5^ and refs. therein) could survive. 
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